Gravitational lensing effects on the gamma-ray burst Hubble diagram 
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Gamma-Ray Bursts (GRBs) offer a potential way to extend the Hubble diagram to very high 
redshifts and to constrain the nature of dark energy in a way complementary to distant type la 
supernovae. However, gravitational lensing systematically brightens distant GRBs through the 
magnification bias, in addition to increasing the dispersions of distance measurements. We investi- 
gate how the magnification bias limits the cosmological usage of GRBs. We perform Monte-Carlo 
simulations of Swift GRBs assuming a cosmological constant dominated universe and then constrain 
the dark energy equation of state neglecting gravitational lens effects. The originally assumed model 
is recovered with 68% confidence limit even when the dispersion of inferred luminosities is compa- 
rable to that of type la supernovae. This implies that the bias is not so drastic for Swift GRBs as 
to change constraints on dark energy and its evolution. However, the precise degree of the bias in 
cosmological parameter determinations depends strongly on the shape of the luminosity function of 
GRBs. Therefore, an accurate determination of the shape of the luminosity function is required to 
remove the effect of gravitational lensing and to obtain an unbiased Hubble diagram. 

PACS numbers: 95.36.+X, 98.70.Rz 
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I. INTRODUCTION 

The expansion history of the universe provides invalu- 
able information about the energy budget. Observations 
of distant type la supernovae (SNela) revealed that the 
universe is accelerating [l|, implying that the significant 
amount of dark energy is present. The nature of dark 
energy is still unknown and is one of central issues in 
modern cosmology. It is expected that detailed measure- 
ments of the expansion history leads to understanding 
the nature of dark energy, because the expansion rate of 
the universe is related with the equation of state of the 
dark energy. 

In addition to SNela, gamma-ray bursts (GRBs) can 
serve as another probe of the expansion history of the 
universe. Although luminosities of different long duration 
GRBs are not quite similar, they can be inferred from 
observables such as variability [3|, spectral-lag Q, peak 
energy |1, Q, and jet opening angle 0- An advantage 
of GRBs is that GRBs are observed at higher redshifts 
than SNela; the average redshift of GRBs discovered by 
Swift satellite is z ~ 3 Q , and GRB with redshift as high 
as z ~ 6.3 was indeed discovered Q- The higher mean 
redshift implies that we may obtain useful information 
on the expansion history from observations of GRBs, in 
a complementary way to SNela |9|. In fact, attempts to 
construct the Hubble diagram from observed GRBs and 
to constrain cosmological parameters have already been 
made jl3|- Given the fact that the number of GRBs is 
now rapidly growing, GRBs are expected to offer unique 
insight into the nature of dark energy. 

However, it should be kept in mind that gravitational 
lensing has sometimes a great impact on high-redshift ob- 
jects. While the effect of gravitational lensing on SNela 
is modest, effectively just introducing additional disper- 
sion (e.g., lU), gravitational lensing is expected to affect 



GRBs much more drastically because of the following two 
reasons. First, the effect of lensing is a strong function of 
the redshift. At higher redshifts, the probability distribu- 
tion functions (PDFs) of lensing magnification has much 
larger dispersions, and also deviates from the Gaussian 
distribution more significantly. Second, and more impor- 
tantly, their apparent isotropic luminosity could range 
over three orders of magnitudes, although they may be 
standardized through the use of luminosity indicators to 
some (uncertain) level. In this situation, the magnifica- 
tion bias, i.e., the effect of observing highly magnified 
GRBs selectively whose fiuxes without lensing effect are 
below a flux limit of GRB detections, becomes quite sig- 
nificant. 

In this paper, we study the effect of gravitational lens- 
ing on the GRB Hubble diagram and the determination 
of the dark energy equation of state. We pay particular 
attention on the non-Gaussian nature of magnification 
PDFs and the magnification bias. First, we derive the 
amount of magnification bias analytically, and then we 
perform Monte-Carlo simulations to show the impact of 
gravitational lensing on the cosmological parameter de- 
termination. As a fiducial cosmological model, we con- 
sider a model with the matter density f^M = 0.27, the 
dark energy density f^DE — 0.73, the dark energy equa- 
tion of state w{z) = —1, the dimensionless Hubble con- 
stant h — 0.72. Below we refer the model as the ACDM. 



II. MAGNIFICATION PROBABILITY 
DISTRIBUTION 

In this section, we derive the magnification PDF which 
is used to examine the effect of gravitational lensing on 
GRBs. We adopt a compound method that we derive the 
PDFs at around the peak and tail separately and combine 
them. First, for the PDF at the peak, we construct it 



from the PDF of lensing convergence k. We adopt a modified log-normal model [13 for the convergence PDF: 
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where Kmin is the convergence of the empty beam, and 
N^ A, Lo are determined from the normalization and the 
conditions J KP{K)dK — and J K^P{K)dK = (k^). Al- 
though the model was intended to describe the PDF for 
each lens plane, we regard it as the PDF for convergence 
projected along all hne-of-sight by adop ting the projected 
variance for (k^)] Taruya et al. [l3] showed that this 
prescription well reproduces the non-Gaussianity of the 
PDF at around the peak, ^/(k^)^/^ < 10. The pro- 
jected variance of lensing convergence, (k^) is computed 
by a standard method using non-linear power spectrum. 
Since GRBs are point sources, we adopt sufficiently small 
smoothing angle 9s = 0.1" in computing the variance. 
We convert it to the magnification PDF by neglecting 
shear and adopt a relation 
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where fj, denotes a magnification factor. On the other 
hand, we compute PDFs at high-magnification tails with 
a halo approach assuming an NFW profile for dark halos 
(see, e.g., 1145). Specifically, we compute the cross sec- 
tion, i.e., the area with the magnification larger than /x, 
a{> fi), assuming an NFW profile [T^, and then they are 
summed up with a mass function of dark halos dn/dM: 



FIG. 1: The magnification PDFs for source redshift z — 1 
(dotted), 3 (solid), and 7 (dashed) are plotted as a function of 
the magnification factor fi. The inset shows expanded view 
of the PDFs around the peaks. 



III. MAGNIFICATION BIAS 
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For the mass function, we adopt a form proposed by 
Sheth & Tormen [I3I . These PDFs are computed assum- 
ing the ACDM model and CTg = 0.8. We connect both 
PDFs where they intersect, which occurs at fj, ^ 2. We 
show the resulting PDFs in Figure^ As seen, the widths 
of the PDFs are wider for higher redshifts. The probabil- 
ities for relatively high magnifications are becoming very 
high at high redshifts. The peaks shift to smaller magnifi- 
cation factor /i with increasing redshift to assure that the 
mean magnification is unity. These behaviors clearly in- 
dicate that effects of gravitational lensing is more promi- 
nent at high redshifts. The PDFs are roughly consistent 
with those in literature (e.g., 0|). We note that the 
detailed shape of the magnification PDF is not very im- 
portant for our purpose; our qualitative result will not 
be affected by a slight change of the PDF. 



In this section, we estimate the impact of the mag- 
nification bias effect analytically. This can be done by 
combining the magnification PDF we derived above and 
GRB luminosity function in the following way. The GRB 
luminosity function (t>(L) has been constrained from the 
number count of GRBs as a function of fluxes and the 
redshift distribution of GRBs H [Hill. We adopt a 
luminosity function derived by Firmami et al. [ig which 
includes the luminosity evolution of GRBs. Specifically 
we adopt a model of a double power-law luminosity func- 
tion with no correlation of peak energy with luminosity (a 
model named 'DPE' in [ia|)- The slopes of the luminosity 
function are —1.53 and —3.4, and the break luminosity 
Lb evolves with redshift such that L^ oc {1 + z)^-^. They 
also constrained the GRB formation rate from the data, 
and pointed out that it resembles the observed cosmic 
star formation rate. 

Given the shape of luminosity function, we can esti- 
mate the effect of the magnification bias. The mean mag- 
nification factor for observed GRBs are calculated from 
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FIG. 2: The mean deviation of the distance modulus (DM) for 
Swift GRBs introduced by gravitational lensing, — 2.51og(^), 
is shown by a thin solid line. The shaded region shows how 
sensitive the magnification bias is to the slope of the lumi- 
nosity function: The width corresponds to the change of the 
slope by ±0.3. Solid and dashed lines are distance moduli in 
evolving dark energy cosmology models with respect to that 
in the ACDM model. For these lines, the gravitational lensing 
effect is not considered. 



the magnification PDF as 
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where Lmin is the minimum luminosity corresponding to 
the flux limit of observations and is a function of the 
redshift. As a specific example, we consider the Swift 
satellite as a GRB detector. Since all GRBs with the 
redshifts measured have observed photon flux larger than 
~ 0.4 cm-^s-i (15 - 150 keV) ^9J, in what follows we 
adopt the flux limit Pnm = 0.4 cm~^s~^. The energy 
band is converted correctly usin g a spectral energy dis- 
tribution of GRBs presented in [ig. The magnification 
factor due to lensing affects the estimate of a distance 
modulus by — 2.51og/i. Thus, we expect that the Hubble 
diagram derived from GRBs is on average shifted from 
the true one by —2.5 log(/i) . Hence this quantity provides 
a good estimate of how much constraints on cosmological 
parameters are biased. 

Figure |2] shows how gravitational lensing biases the 
Hubble diagram from GRBs. First of all, the degree of 
the bias again depends strongly on the redshift. At z = 5 
lensing decreases the distance modulus by ~ —0.05, thus 



it cannot be ignored. Second, errors of slopes of luminos- 
ity functions introduce large uncertainties in the effect 
of gravitational lensing on the GRB Hubble diagram, 
in particular at high redshifts. The shallow luminosity 
function means that there are much more GRBs below 
the flux limit, thus we have much more GRBs that are 
beyond the flux limit just because of ampliflcations by 
gravitational lensing than those which are below the flux 
limit due to dimming by lensing, resulting in a large mean 
magnification. The dispersion introduced by lensing is 
^ 0.3 mag at z ~ 5, larger than the systematic shift by 
the magnification bias. In summary, we find that gravi- 
tational lensing systematically changes the shape of the 
GRB Hubble diagram particularly at high redshifts, not 
to mention introducing additional dispersions. 

How does the bias affect the determination of the dark 
energy equation of state from the GRB Hubble diagram? 
To explore this, we consider the following two parame- 
terizations of the equation of state w{z): One is w{z) = 
wq -\- [z/(l -t- z)]wa and the other is w{z) — wq + w' z, 
both of which have been widely adopted in dark energy 
studies. We note that in the latter parametrization w{z) 
increases rapidly with increasing redshift, which results in 
unphysical values oi w{z) at high redshifts. We compute 
distance moduli for following two parameter sets, (r^DE, 
Wo, w')=i0.75, -1.1, 0.4) and {^be, wq, Wa)={0.74, 
— 1.2, 1.0), which are plotted in Figure 13 As shown in 
the Figure, distance moduli of these evolving dark energy 
models behave just as that of the ACDM model biased 
by gravitational lensing. The increase of the equation 
of state as a function of z results in the enhancement 
of the expansion rate at high-redshifts and thus make 
distant objects brighter than those in no dark energy 
evolution model, which is just similar as gravitational 
lensing effect. It is interesting to note that the behavior 
more resembles dark energy parametrized unphysically 
as w{z) — Wo +w'z. Therefore the effect of gravitational 
lensing is expected to appear as an artificial evolving dark 
energy in determining cosmological parameters from the 
GRB Hubble diagram. 



IV. MONTE-CARLO SIMULATIONS 

We pursue the bias introduced by gravitational lensing 
further by performing Monte-Carlo simulations of GRB 
observations. Again, we adopt a model of Firmami ct 
al. [I3I for the cosmological distributions of GRBs, and 
also assume the ACDM model as the underlying cos- 
mology. First we distribute GRBs in redshift-luminosity 
space according to the GRB formation rates and the best- 
fit luminosity function. We are conservative to restrict 
our attention to GRBs at z < 7, because it is still un- 
clear whether we can really measure redshifts of GRBs 
at z > 7. For each event, we first compute an observed 
flux. After adding a measurement error which is modeled 
by a Gaussian error with a = 0.05 mag, we compare it 



with the flux limit of Swift satellite Pn 



0.4 cm ^s 




FIG. 3; Projected constraints on the dark energy equation of 
state expected from 400 Swift GRBs. Light and dark shaded 
regions indicate 68% limits for the cases with luminosity in- 
ference uncertainties of cr = 0.6 mag and 0.16 mag, respec- 
tively. The best-fit values for a = 0.6 is shown by crosses. 
Filled square indicates the ACDM model that we originally 
assumed in simulations of GRBs. For the parameterization of 
the equation of sate, we adopt w{z) = wo + [z/i^ + z)]'Wa- 



and select it if the flux is above the limit. In this way 
we construct simulated data of 400 GRBs; the number is 
based on the fact that Swift satellite is discovering ~ 100 
GRBs per year. We note that the number is somewhat 
optimistic since only of order half of Swift GRBs have 
measured redshifts. Next, for each event, we infer its 
absolute luminosity using any correlations with observ- 
ables such as variability and spectral-lag. Here we do not 
ask which indicators are used, but just assume a Gaus- 
sian with a = 0.6 mag as an error associated with the 
estimate of the absolute luminosity. This error is conser- 
vative in the sense that one can reduce the error to ^ 0.45 
mag by making use of several indicators (see [23|), but 
could be optimistic in the sense that it ignores several 
systematic effects [23|. It is possible that we find better 
correlations which significantly reduce the error in esti- 
mating the absolute luminosity. Thus, we also consider 
the case a = 0.16 mag which is comparable to that of 
SNela. From observed fluxes and inferred absolute lumi- 
nosities, we can derive distance moduli for these GRBs, 
which are used to constrain cosmological parameters. We 
compute x^ in four parameter space, Om = 1 — f^DE, h, 
wq, and Wa assuming a flat universe. We add Gaussian 
priors ^m = 0.27 ± 0.04 and h = 0.72 ± 0.04, which are 
reasonable given the current accuracy of cosmological pa- 
rameter determinations |22] . To reduce offsets of best- fit 
values by statistical fiuctuations, we repeat the calcula- 
tion of x^ foi' 30 different realizations of 400 GRBs and 
take an average of x^ ■ 



We show constraints on the equation of state in Figure 
13 As expected, gravitational lensing shifts the best-flt 
models: The best fit parameter set from fitting is (wq, 
Wa)=(— 1.15, 0.6). However, the assumed model is still 
within 68% confidence limit for both large and small er- 
rors in estimating GRB absolute magnitudes, implying 
that gravitational lensing does not bias the parameter 
estimate very much for Swift GRBs. The cosmological 
constant model (wiz) = — 1 has Ax^ ~ 0.5 when the 
error is large and Ax^ ~ 0.1 when the error is small. 
We note that the effect of lensing is dependent on the 
shape of the luminosity function: Thus the effect could 
be much more significant if the luminosity function has a 
steep slope at high redshifts, as is clear from Figure|21 In 
summary, the gravitational lensing is important system- 
atic effect in constraining the nature of dark energy from 
the GRB Hubble diagram but is probably not so drastic 
as to totally change results from Swift GRBs. 



V. CONCLUSION 

In this paper, we have studied the effect of gravita- 
tional lensing on the determination of dark energy prop- 
erties from distant GRBs. Gravitational lensing sys- 
tematically brighten the apparent luminosity of observed 
GRBs, resulting in the modification of the Hubble dia- 
gram inferred from GRBs. Our Monte-Carlo simulations 
have shown that lensing does bias the best-fit model, but 
we can recover the assumed model within 68% confidence 
limit for Swift GRBs. Therefore, it does not degrade the 
cosmological usage of GRBs very much, at least within 
accuracies currently expected from Swift observations. 
One of the reasons is that the constraint on dark en- 
ergy from Swift GRBs are not very strong. We, how- 
ever, emphasize that the amount of the bias introduced 
by lensing is quite sensitive to the shape of the luminos- 
ity function. Thus, we need an accurate measurement 
of the GRB luminosity function, in particular near the 
flux limit and beyond, in order to remove the effect of 
gravitational lensing and to obtain unbiased Hubble di- 
agram. In addition, the strong influence of gravitational 
lensing on high-redshift GRBs implies that estimates of 
high-redshift GRB (or similarly star formation) rates and 
even the luminosity function itself may be biased. We 
plan to address these issues in a forthcoming paper. 
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